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MMSE  EQUALIZATION  FOR  AERONAUTICAL 
TELEMETRY  CHANNELS 


Michael  Rice 

Brigham  Young  University 

Md.  Shah  Afran 

Mohammad  Saquib 

The  University  of  Texas  at  Dallas 


ABSTRACT 

This  paper  presents  performance  analysis  of  the  minimum  mean  squared  error  (MMSE)  equalizers 
applied  to  aeronautical  telemetry  channels.  The  challenge  for  equalizing  received  samples  of  the 
modulated  signal  lies  in  the  fact  that  the  underlying  continuous-time  SOQPSK-TG  waveform  is 
not  wide-sense  stationary.  However  it  is  assumed  so  in  order  to  meet  real-time  implementation 
requirements.  Two  approximations  of  the  autocorrelation  function  of  the  SOQPSK-TG  waveform 
are  used  for  designing  MMSE  equalizers.  Their  performance  are  investigated  against  the  zero 
forcing  equalizer  for  measured  aeronautical  telemetry  channels. 


INTRODUCTION 

The  propagation  of  the  radio  signal  from  an  airborne  transmitter  to  a  ground-based  receiver  over 
multiple  paths  may  cause  multipath  interference.  Usually,  one  of  the  paths  is  the  line-of-sight 
propagation  path  whereas  the  others  are  due  to  reflections.  Multipath  interference  continues  to  be 
the  dominant  cause  of  link  outages  in  aeronautical  telemetry.  In  this  paper  we  investigate  a  data- 
aided  approach  to  equalization  assuming  iNET  packet  structure.  In  data-aided  equalization,  the 
equalizer  filter  coefficients  may  be  computed  from  the  multipath  channel  coefficients. 

iNET-formatted  transmissions  include  a  128-bit  preamble  and  64-bit  attached  sync  marker  (ASM) 
preceding  a  block  of  data  bits  (at  least  6144  bits:  an  LDPC  codeword):  see  Figure  1.  Since 
the  preamble  and  ASM  bits  are  known,  the  receiver  can  compare  the  received  signal  to  a  locally 
stored  copy  of  the  SOQPSK-TG  signal  corresponding  to  the  preamble  and  ASM  bit  fields.  This 
comparison  is  capable  of  producing  estimates  of  the  frequency  offset,  noise  variance,  and  multipath 
channel  coefficients  [4].  The  multipath  channel  coefficient  estimates  can  then  be  used  to  obtain 
equalizer  filter  coefficients. 
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PRE 

ASM 

DATA 

(128  bits) 

(64  bits) 

(6144  bits) 

Figure  1:  The  iNET  packet  structure  used  in  this  paper. 


The  minimum  mean-squared  error  (MMSE)  filter  coefficients  depend  on  multipath  channel  coeffi¬ 
cients,  autocorrelation  function  of  the  SOQPSK-TG  waveform  and  noise  variance.  Unfortunately, 
SOQPSK-TG  waveform  is  not  wide-sense  stationary  which  results  in  time-varying  filter  coeffi¬ 
cients  and  makes  MMSE  equalizer  practically  very  difficult  to  realize.  This  situation  leads  us 
to  make  two  approximations  of  the  autocorrelation  function  of  the  SOQPSK-TG  waveform  and 
investigate  their  performance  against  the  zero  forcing  (ZF)  equalizer  for  measured  aeronautical 
telemetry  channels. 


SYSTEM-LEVEL  DESCRIPTION 


The  bit  sequence  for  iNET  is  depicted  in  Figure  1.  The  preamble  sequence  (PRE)  is  CD98hex  re¬ 
peated  eight  times  [5,  p.  48].  The  preamble  field  is  followed  by  the  attached  sync  marker  (ASM) 
field  defined  as  034776C7272895B0hex.  The  DATA  field  is  6144  randomized  data  bits.1  The  iNET 
bit  sequence  is  modulated  by  SOQPSK-TG  waveform  which  propagates  through  a  frequency  selec¬ 
tive  channel  and  experiences  a  frequency  offset  as  well  as  the  addition  of  additive  white  Gaussian 
noise. 


The  received  signal  is  filtered,  down-converted  to  I/Q  baseband,  and  sampled  (not  necessarily  in 
that  order)  using  standard  techniques.  The  resulting  sequence  of  received  samples  is 


r{n ) 


'  n2 

^  h{k)s{n  —  k ) 

_fc=-AT1 


ej“on  +  w{n), 


(1) 


where  h(n )  is  the  impulse  response  of  the  equivalent  discrete-time  channel  with  support  on  —N{  < 
n  <  N2,  w0  rads/sample  is  the  frequency  offset,  and  w(n)  is  a  complex- valued  zero-mean  Gaussian 
random  process  with  variance  rU  . 

The  focus  of  this  paper  is  on  equalizing  the  I/Q  baseband  samples  of  the  received  signal  (1).  Prior 
to  applying  the  equalization  techniques  multiple  tasks  need  to  be  performed  by  the  receiver  as 
shown  in  Figure  2.  The  preamble  and  ASM  bits  are  known  and  thus  the  samples  corresponding  to 
those  bits  are  used  to  estimate  the  frequency  offset,  channel  impulse  response,  and,  for  the  MMSE 
equalizer,  the  noise  variance.  Before  these  estimations  can  be  performed,  the  start  of  the  samples 
corresponding  to  the  preamble  bits  in  the  received  signal  must  be  detected.  This  is  accomplished 
by  the  preamble  detector  block,  whose  algorithm  is  based  on  the  detection  algorithm  described  in 
[7].  Once  the  start  of  the  preamble  is  known,  the  frequency  offset  is  estimated  using  the  algorithms 
described  in  [4].  The  frequency  offset  is  used  with  a  complex-exponential  to  derotate  the  received 

'These  bits  correspond  to  a  single  LDPC  codeword  in  the  coded  system.  Here,  we  evaluate  the  uncoded  bit  error 
rate  (BER)  after  equalization. 
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Figure  2:  The  data  packet  format  and  high-level  signal  processing  explored  in  this  paper. 


data  to  remove  the  frequency  offset.  The  derotated  data  rd(n )  are  used  to  estimate  the  channel  and 
noise  variance  as  described  in  [4].  The  channel  estimates  h(n),  for  —N\  <  n  <  N2,  are  then  used 
to  compute  the  MMSE  and  ZF  equalizer  filter  coefficients. 


THE  EQUALIZATION  ALGORITHMS 

Since  SOQPSK-TG  is  a  nonlinear  modulation,  the  equalizer  cannot  operate  on  the  symbols  in  the 
same  way  it  does  for  linear  modulation  (cf.,  [6,  Chapter  9]).  Consequently,  the  equalizer  must 
operate  on  the  samples  of  SOQPSK-TG,  similar  to  the  way  fractionally  spaced  equalizers  operate. 
The  equalizers  operate  in  the  system  configuration  shown  in  Figure  3  [cf.,  Figure  2].  Here,  the 
derotated  samples  r<i{n )  are  equalized  using  a  length  L\  +  L2  +  1  FIR  filter  defined  by  the  impulse 
response  c(n )  for  —  Li  <  n  <  L2  to  produce  the  output 

l2 

y(n)  =  c(m)rd(n  —  m).  (2) 

m=—Li 

The  equalizer  output  forms  the  input  to  the  well-known  symbol-by-symbol  SOQPSK  detector  com¬ 
prising  a  detection  filter  operating  at  N  =  Tb/T  samples/bit  and  a  decision  process,  operating  on 
the  decision  variable  u(k )  at  1  sample/bit.  This  detector,  based  on  an  offset  QPSK  approximation 
of  SOQPSK-TG,  is  described  in  more  detail  in  [8,  9].  The  detectors  of  Figure  3  also  include  a 
phase  lock  loop  (PLL).  The  PLL  is  required  to  track  out  any  residual  phase  increments  due  to 
frequency  offset  estimation  errors.  A  timing  loop  is  not  required  because  timing  offsets  are  part 
of  the  channel  estimate  h{n).  Now  we  will  organize  the  MMSE  equalizer  filter  coefficients  into 
(L1  +  L2  +  1)  x  1  vectors  as  follows: 


cmmse(~  L\ 


CMMSE  — 


Cmmse(O) 


(3) 


L  cMMSe(-^2)  J 
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The  MMSE  equalizer  is  a  filter  that  minimizes  the  mean  squared  error 


S 


E 


/  s(n)  —  rain)  *  c{n) 


(4) 


As  mentioned  earlier,  the  challenge  with  equalizing  samples  of  the  modulated  signal  is  that  the 
underlying  continuous-time  waveform  is  not  wide-sense  stationary  [6].  This  fact  carries  over  by 
the  autocorrelation  function  of  s(n ) 

Rs(k,£)  =  ^E{s(k)s*(£)J.  (5) 

Notice  in  (5)  that  the  autocorrelation  function  is  a  function  of  both  sample  indexes,  not  the  dif¬ 
ference  between  them.  Consequently,  the  equalizer  filter  coefficients  are  a  function  of  the  sample 
index  n.  It  is  hard  to  see  how  this  solution  has  any  practical  utility,  especially  in  the  presence  of  a 
real-time  performance  requirement.  In  the  end,  the  designer  is  left  with  suboptimal  approaches  of 
reduced  computational  complexity  whose  accompanying  performance  penalty  is  acceptable. 

The  simplest  suboptimal  approach  is  to  assume  that  the  signal  samples  are  wide-sense  stationary. 
Here,  the  autocorrelation  function  is  of  the  form 

Rs(k-£)  =  ^E{s(k)s*  (£)};,  (6) 

that  is,  the  autocorrelation  function  depends  on  the  difference  of  the  sample  time  indexes.  The 
wide-sense  stationary  assumption  for  s(n)  greatly  simplifies  the  solution.  Because  the  equalizer 
coefficients  no  longer  depend  on  the  samples  index  n,  the  relationship  between  s(n)  and  the  equal¬ 
izer  output  s(n )  is 

L2 

s(n )  =  c{n)  *  r<i{n )  =  c{m)rd{n  —  m) .  (7) 

m=—L\ 

Recall  that  rd{n )  is  the  derotated  version  of  the  received  samples.  The  vector  of  filter  coefficients 
that  minimizes  the  mean  squared  error 


S 


E 


(8) 


is  given  by 


c  = 


GRc  iGf  +  Ru 


Rs,2gf 


(9) 


where  c  is  the  (L 1  +  L2  +  1)  x  1  vector  of  filter  coefficients,  G  is  the  ( L\  +  L2  +  1)  x  (N\  +  N2  + 
Li  +  L2  +  1)  matrix  described  by 


h(N2)  •••  hi-Ni) 

h(N2)  •••  h(~N 1) 

h(N2)  •••  h(-Ni). 


(10) 
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Rs,i  is  the  (Li  +  L2  +  Ari  +  A^  +  1)  x  (A  1  +  A2  +  Ari  +  N2  +  1)  matrix 


Rs,i  — 


Rs(  0) 

Rs(l) 


Rs(-l) 
Rs(  0) 


Rsi-Lt  -  L2  -  m  -  N2) 
R,s(—L\  -L2-  N1^N2  +  1) 


RS{L  1  +  A2  +  A^i  +  A^2)  Rs(Li  +  A2  +  A'i  +  A^2  —  1)  •  •  • 

R,  is  the  (Li  +  L2  +  1)  x  (Lx  +  L2  +  1)  noise  autocorrelation  matrix  given  by 

-Ro(O)  •••  Rw(—Li  —  L2) 

R  ;  ; 

Rw(Li  +  A2)  •  •  •  -Rui(O) 

Rs  2  is  the  ( Li  +  L2  +  1)  x  ( Li  +  L2  +  1)  matrix  given  by 


R 


■s,  2  — 


Rs(0) 

Rs(l) 


Rs(- 1) 
Rs(0) 


Rs(Ai  +  L2)  Rs(Li  +  L2  —  1) 


Rs{— L\  —  L2) 
Rs(—L\  —  L2  +  1) 

Rs(  0) 


Ra{  0) 


(ID 


(12) 


(13) 


and  g  is  the  1  x  ( L1  +  A2  +  1)  vector  given  by 

g  =\h(Li)  •"  h(-L2)\,  (14) 

where  it  is  understood  that  h(n)  =  0  for  n  <  —N\  or  n  >  N2  (how  many  zeros  need  to  be 
prepended  and  appended  depends  on  the  relationship  between  Li  and  N2  and  the  relationship 
between  L2  and  A^). 

The  question  is  now,  what  function  should  be  used  for  the  autocorrelation  function  Rs(k)l  Two 
approximations  are  investigated  here.  The  first  is  an  empirically-derived  autocorrelation  function. 
The  empirical  autocorrelation  function  is  obtained  by  generating  a  large  number  of  samples  s(n) 
and  using  the  standard  estimation  technique  assuming  wide  sense  stationarity.  Given  L  samples  of 
s(n )  for  n  —  0, 1, .. . ,  L  —  1,  this  empirical  autocorrelation  function  is 


Re{k) 


1 

2  (L-k) 


L- 1 

^2  s(n)s*(n  —  k ), 

n=k 


0  <  k  <  L-  1 


(15) 


together  with 


Re(k)  =  R*(—k),  —L  <  k  <  0. 


(16) 


A  plot  of  Re(k )  corresponding  to  L  =  2  x  106  samples  of  SOQPSK-TG  sampled  at  2  samples/bit 
is  shown  in  Figure  4  for  the  first  100  lags  (i.e.,  —100  <  k  <  100).  The  top  plot  shows  the  real 
part  of  Re(k)  and  the  lower  plot  shows  the  imaginary  part  of  Re(k).  The  only  significant  values 
are  those  for  —5<k<5  and  indicated  by  markers  on  the  plot.  Consequently,  in  the  simulation 
results  presented  below,  we  assume  Re(k )  =  0  for  \k\  >5. 
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Figure  4:  A  plot  of  the  empirical  autocorrelation  function  for  SOQPSK-TG:  (top)  the  real  part  of 
Re(k)-  (bottom)  the  imaginary  part  of  Re(k).  The  sample  rate  for  the  SOQPSK-TG  samples  is  at 
2  samples/bit.  Markers  indicate  the  values  for  —  5  <  k  <  5. 
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frequency  Gaussian 
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Figure  5:  A  block  diagram  of  the  simulation  procedure. 


The  second  approximation  is  to  assume  the  data  are  uncorrelated, 
function  of  the  form 


Ri(k)  =  a2sS(k). 


This  generates  a  correlation 


(17) 


Here,  the  corresponding  correlation  matrices  Rs q  and  Rs  2  are  scaled  identity  matrices  and  they 
function  as  regularizes  in  the  numerical  computations  (9). The  solution  is  given  by  a  form  of  the 
Wiener-Hopf  equations  [10].  The  +  L2  +  1)  x  1  vector  of  MMSE  equalizer  filter  coefficients 
are 


Cmmse  — 


GGf  +  ^flll+l5+1 

cri 


-l 


(18) 


where  G  is  the  {L\  +  L2  +  1)  x  {N\  +  N2  +  Li  +  L2  +  1)  matrix  as  described  in  (10),  g  is  the 
1  x  (Li  +  Lo  +  1)  vector  as  described  in  (14)  and  a2  is  the  signal  power.  This  solution  is  based 
on  the  assumption  that  SOQPSK-TG  samples  corresponding  to  a  sample  rate  of  2  samples/bit  are 
approximately  uncorrelated.  The  ZF  equalizer  is  a  filter  that  is  the  best  length- (Li  +  L2  +  1)  FIR 
approximation  to  the  “inverse”  of  the  channel.  Its  filter  coefficients  can  be  derived  from  (18)  by 
letting  noise  variance  b2w  =  0. 


PERFORMANCE  RESULTS 

The  BER  performance  of  the  equalization  techniques  was  assessed  using  the  simulation  environ¬ 
ment  outlined  in  Figure  5.  The  simulation  parameters  were  the  following: 

1 .  The  payload  data  rate  was  equivalent  to  10  Mbits/s  (the  equivalent  “over-the-air”  bit  rate  was 
10.3125  Mbits/s).  The  iNET-formatted  SOQPSK-TG  signal  and  channel  were  generated  at 
an  equivalent  sample  rate  of  2  samples/bit. 

2.  Because  the  channel  estimator  does  not  know  the  true  length  of  the  channel,  the  estimator 
used  values  for  N\  and  N2  larger  than  any  of  the  test  channels.  These  values  were  Ar,  =  12 
and  N2  =  25  samples. 

3.  The  equalizers  used  L\  =  4  x  N\  =  48  samples  and  L2  =  4  x  jV2  =  100  samples.  Thus  the 
length  of  equalizer  filter  was  149  samples. 
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Table  1:  Description  of  the  ten  test  channels  used  in  the  simulations. 


channel 

N, 

n2 

length 

environment 

01 

1 

7 

9 

Taxiway  E 

02 

2 

17 

20 

Taxiway  E 

03 

1 

22 

24 

Taxiway  E 

04 

5 

13 

19 

Takeoff  on  22L 

05 

1 

1 

3 

Cords  Road 

06 

1 

2 

4 

Cords  Road 

07 

0 

4 

5 

Cords  Road 

08 

2 

3 

6 

Black  Mountain 

09 

1 

1 

3 

Black  Mountain 

10 

2 

3 

6 

Land  on  22L 

4.  The  simulations  were  performed  over  10  representative  channels  derived  from  channel  sound¬ 
ing  measurements  conducted  at  Edwards  AFB  under  the  M4A  program  [12].  The  test  chan¬ 
nels  are  summarized  in  Table  1  and  the  corresponding  frequency-domain  plots  are  shown  in 
Figure  6. 

The  simulated  BER  performance  is  shown  in  Figures  7  -  1 1 .  In  all  cases  we  observe  that  both  ver¬ 
sions  of  the  MMSE  equalizer  exhibit  almost  identical  performance,  which  is  noticeably  better  than 
the  performance  of  the  ZF  for  most  of  the  channels  except  channels  6  and  9.  This  is  to  be  expected 
because  test  channels  6  and  9  are  rather  benign  [see  Figure  6].  Similarly,  when  the  channel  has  deep 
and  wide  spectral  nulls  in  the  middle  of  the  spectrum  of  SOQPSK-TG  waveform,  MMSE  equalizer 
is  expected  to  provide  significant  gain  over  the  ZF  equalizer.  In  Figure  8  for  channel  4,  it  can  be  no¬ 
ticed  that  at  a  target  BER  =  10-5  both  MMSE  equalizers  yield  about  3  dB  signal-to-noise  ratio  gain 
over  the  ZF  equalizer.  The  fact  that  the  ZF  equalizer  has  such  performance  is  not  surprising.  The 
ZF  equalizer  simply  “inverts  the  channel”  (this  phrase  derives  from  the  frequency  domain  point 
of  view).  For  channels  with  nulls,  the  “inversion”  restores  the  frequency  content  of  the  desired 
signal  in  the  frequency  band  surrounding  the  null.  This  restoration  also  amplifies  the  noise  power 
in  the  same  frequency  band.  The  end  result  is  a  phenomenon  known  as  “noise  amplification:”  the 
distortion  due  to  the  multipath  channel  is  corrected  at  the  cost  of  reduced  signal-to-noise  ratio.  In 
contrast  an  MMSE  equalizer  takes  a  more  measured  approach  to  “channel  inversion”  and  balances 
the  impact  of  residual  multipath  distortion  and  amplified  noise  on  cost  functions  [mean  squared 
error  (4)  for  the  MMSE  equalizer].  Again,  from  the  plots  of  simulated  BER  we  observe  that  there 
is  essentially  no  difference  in  the  performance  between  the  two  versions  of  the  MMSE  equalizer. 
Consequently,  R{(k)  is  preferable  over  Re(k)  because  this  choice  simplifies  the  computations  of 
the  equalizer  filter  coefficients. 
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Figure  6:  Frequency-domain  plots  of  the  example  channels  from  channel  sounding  experiments  at 
Edwards  AFB.  In  each  plot,  the  thick  line  is  the  channel  frequency  response  and  the  thin  line  is  the 
power  spectral  density  of  SOQPSK-TG  operating  at  10.3125  Mbits/s. 
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Figure  7:  Simulation  results  for  test  channel  1  (left)  and  test  channel  2  (right). 


Figure  8:  Simulation  results  for  test  channel  3  (left)  and  test  channel  4  (right). 
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Figure  9:  Simulation  results  for  test  channel  5  (left)  and  test  channel  6  (right). 


Figure  10:  Simulation  results  for  test  channel  7  (left)  and  test  channel  8  (right). 
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Figure  11:  Simulation  results  for  test  channel  9  (left)  and  test  channel  10  (right). 

CONCLUSIONS 

This  paper  demonstrated  the  effectiveness  of  the  MMSE  equalizers  against  the  ZF  eqalizer  with 
iNET-formatted  SOQPSK-TG  over  measured  aeronautical  telemetry  channels.  SOQPSK-TG  wave¬ 
form  is  not  wide-sense  stationary.  However,  the  real-time  implementation  of  the  MMSE  equalizer 
led  us  to  assume  the  underlying  continuous  waveform  as  a  wide-sense  stationary  process.  To 
further  simplify  the  computational  complexity  of  the  equalizer  coefficients  we  assumed  that  the 
waveform  was  not  only  wide-sense  stationary  but  also  uncorrelated.  Our  numerical  results  showed 
that  this  simplification  did  not  penalize  the  system  in  terms  of  BER  and  depending  on  the  channel, 
at  a  target  BER  =  10-5  MMSE  equalizers  were  capable  of  providing  0  —  3  dB  signal-to-noise  ratio 
gain  over  the  ZF  equalizer. 
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